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Modified Associate Formalism without Entropy Paradox: 
Part I. Model Description 
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Abstract 

A Modified Associate Formalism is proposed for thermodynamic modelling of solution phases. The approach is free from the entropy paradox 
described by Luck et al. (Z. Metallkd. 80 (1989) pp. 270-275). The model is considered in its general form for an arbitrary number of solution 
components and an arbitrary size of associates. Asymptotic behaviour of chemical activities of solution components in binary dilute solutions 
is also investigated. 
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1. Introduction 

' The associate model in its various modifications has been 
successfully used for modelling solution phases of metallur- 
gical and chemical engineering interest [1-7]. Besmann and 
Spear [3], for example, utilised the modified associated species 
model for glasses used in nuclear waste disposal. Recently, 
Yazhenskikh et al. [5-7] have successfully applied the asso- 
ciate species model to model melting behaviour of coal ashes 
which is an important problem in coal gasification technolo- 
gies. Good agreement between model predictions and available 
experimental data was reported. 

According to the classical associate model described by Pri- 
gogine and Defay [8], the strong interactions result in the for- 
mation of stable configurations of mixing particles, the so- 
called association complexes or, briefly, associates. Those par- 
ticles that are not involved in the formation of associates are 
called free particles or, interchangeably, monoparticles. The 
associated solution is then considered to be an ideal solution 
of monoparticles and different associates. For example, a bi- 
nary associated solution of components A and B, in which 
only AB-associates are formed, is considered to be a ternary 
ideal solution of the A-monoparticles, B-monoparticles and AB- 
associates. 



The Gibbs free energy G of the associated solution of n A 
moles of the solution component A and n B moles of B is then 
given by 



G = n Al g Al + n B] g Bl +n AB g AB - TS con f, 



(1) 



where the configurational entropy of mixing S con f is expressed 
as 

■Sconf = R i n A \ ln*A, + n Bl \nx Bl + n AB lnx AB ) . (2) 

Here, n A[ , n B[ , n AB are the mole numbers and g A[ , g B[ , g AB 
are the molar Gibbs free energies of the A-monoparticles, B- 
monoparticles and AB-associates, respectively; T is the abso- 
lute temperature and R is the universal gas constant. The molar 
fractions x Al , x B[ and x AB are defined in a usual way. For ex- 
ample, x A] = n Al I (n A[ +«B[ +n AB ). The other molar fractions 
are defined similarly. 

The equilibrium values of the mole numbers n A[ ,n B{ and n AB 
are determined by minimising the Gibbs free energy G given 
by Eq. (1), subject to the mass balance constraints 



n A = n Al +n AB , 
n B = n Bl +n AB . 



(3) 



The adjustable parameter of the model is the molar Gibbs free 
energy Ag AB of the reaction 



A+B^AB. 



(4) 
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Besmann and Spear [3] considered an ideal mixture of as- 
sociate species (instead of monoparticles and associates). The 
stoichiometry of the associate species was specified so that all 
the species contain two non-oxygen atoms per formula unit. In 
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Fig. 1. The configurational entropy of mixing for the associated solution, 
in which only the AB-associates are formed, in the high temperature limit 
compared with the ideal entropy of mixing. 

this approach, contributions of different species to the config- 
urational entropy of mixing are equally weighted. 

Luck et al. [9] described a remarkable feature of the config- 
urational entropy of mixing given by Eq. (2). This feature was 
referred to as an entropy paradox and is briefly described in the 
next section. 

2. Entropy paradox 

Luck et al. considered the high temperature limit for the as- 
sociated solution. The temperature T is assumed to be so high 
that the entropy term plays a dominant role in the Gibbs free 
energy of the solution and the enthalpy changes on forming 
different associates can be neglected. The configurational en- 
tropy of mixing given by Eq. (2) in the high temperature limit 
is higher than that used in the regular solution model (the en- 
tropy of ideal mixing of the solution components). For exam- 
ple, the configurational entropy of mixing of the associated so- 
lution in which only the AB-associates are formed is compared 
with the ideal entropy of mixing in Fig. 1. At the same time, 
the formation of associates models short range ordering in the 
solution. As pointed out by Luck et al. [9], it is a paradoxical 
result that the configurational entropy, which is a measure of 
disorder, appears to be higher in a solution with ordering than 
in a completely disordered solution. 

Another interpretation of this entropy paradox was described 
by Pelton et al. [10]. Consider the binary associated solution 
where only AB-associates are formed and assume that there are 
no Gibbs free energy changes on forming AB-associates from 
monoparticles, so that AgAB = 0. In this case, the configura- 
tional entropy of mixing of the solution should be equal to that 
of an ideal solution, since no interactions between mixing par- 
ticles are assumed. Eq. (2) however, leads to higher values for 
the configurational entropy of mixing that reduces to the ideal 
configurational entropy of mixing only when AgAB = +°°; see 
also Fig. 1. The overestimation of the configurational entropy 



can result in either underestimation of the non-configurational 
entropy or overestimation of the enthalpy of mixing or both. 
This, in tern, can undermine the predictive capabilities of the 
model. 

As pointed out by Luck et al. [9] and later by Pelton et 
al. [10], the expression for the configurational entropy of mix- 
ing used in the quasichemical model does reduce to the ideal 
configurational entropy of mixing when AgAB = 0. In the next 
section, we propose the Modified Associate Formalism which 
is free of the entropy paradox. The paradox is resolved by dis- 
tinguishing between all possible spatial arrangements of parti- 
cles in an associate that have not been taken into account in 
previous associate models. 

3. Model assumptions 

The model proposed in this paper is based on the following 
assumptions. 

1) Similarly to the classical associate model [8], we assume 
that interactions between mixing particles result in the 
formation of associates which are in a stable dynamic 
equilibrium with each other. The associates of the model 
are understood as a tool for modelling short-range inter- 
actions between mixing particles. The associates of the 
model, however, may represent real associated complexes 
present in solution phases. 

2) We assume that the associates do not interact with each 
other and are uniformly distributed (ideally mixed) over 
a lattice. Equivalently, the occupancies of the sites of the 
associate lattice are stochastically independent and have 
identical probability distributions. 

3) In contrast to the classical model, all pure solution com- 
ponents and the chemical solution of these components 
are treated in a unified way. More precisely, we assume 
that the solution and all its components consist of nonin- 
teracting associates of the same size, so that the associates 
are composed of the same number of particles. In this 
approach, the contributions of different associates to the 
configurational entropy of mixing are equally weighted. 

4) Following the convention (see, for example, Ref. [11] 
for more details), we also assume that particles of the 
same type are indistinguishable, while particles of differ- 
ent types and particle sites within an associate are distin- 
guishable. 

4. Model description 

For simplicity of exposition, we exemplify the model by 
considering a binary solution A — B where the associates are 
composed of three particles. Let iia and ris be the mole numbers 
of A and B particles in the solution. Since, by assumption 4) 
above, the particle sites within an associate are distinguishable, 
we also assume that they are numbered. If the 1st and 2nd sites 
in an associate are both occupied by A-particles, while the 3rd 
one is occupied by a B-particle, such an associate is said to be 
of type [AAB] . Other types of associates are defined similarly. 
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Thus, in the solution considered, there are 2 3 — 8 different types 
of associates, [AAA], [BAA], [ABA], [AAB], [BBA], [BAB], [ABB], 
[BBB]. It is important to note that all these types should be 
taken into account in calculating the configurational entropy. 

Let be the mole number of [ijk] -associates, where the 
triplet of symbolic indices i,j,k —A,B specify the associate 
type. Then the mass balance constraints take the form 



n A = 3«[AAA] 

+ 2 (« [BAA] + n [ABA] + n [AAB] ) 
+ n [BBA] + n [BAB] + n [ABB] , 

n B = 3«[ BBB j 

+ 2 (« [ABB] + n [BAB] + n [BBA] ) 
+ n [AAB] + n [ABA] + n [BAA] ■ 



(5) 



By a standard combinatorial argument, the number of available 
microstates £2 is 



£2 



No I n [ijk] ! 

\ i,j.k=A,B UI J 

n (Non m )\ ' 

i,j,k=A,B 



(6) 



where N a is the Avogadro number. Under assumption 2) of 
ideal mixing of associates, the configurational entropy of mix- 
ing S conf is 

Sconf = k B \nQ. , (7) 
where Icr is Boltzmann's constant. Applying the Stirling for- 
mula to Eqs. (6) and (7), 



S C onf=-fi E n [ijk]^ nx [ijk]i 
i,j,k=A,B 



(8) 



where 



*[ijk] - 



L n \i'fk'] 

i',f,k'=A,B 



is the molar fraction of the [ijk] -associates. Therefore, the Gibbs 
free energy of the solution is given by 



G = E n m](8[ijk]+RT\nx m ) 

i,j,k=A,B 



(9) 



where g^] is the molar Gibbs free energy of [ijk] -associates. 
In one mole of pure solution component A there are 1/3 moles 
of [AA4]-associates. Hence, g[AAA] = 3gA, where gA is the molar 
Gibbs free energy of the pure solution component A. Similarly, 
8 [bbb] = 3gB, where g B is the molar Gibbs free energy of the 
pure solution component B. 

The associates, which consist of both A and B particles, will 
be referred to as mixed associates. For the binary solution con- 
sidered, there are 8—2 = 6 types of mixed associates. Their 
molar Gibbs free energies g[ijk], with [ijk] ^ [AAA], [BBB], are 
adjustable parameters of the model. These parameters can de- 
pend on the temperature T as 



g[ijk] = h[ ijk ] 



Ts m , 



(10) 



where fyy^j and are the molar enthalpy and entropy of [ijk]- 
associates, respectively. Alternatively, as adjustable parameters 



of the model, one can employ the Gibbs free energies of the 
reactions of forming mixed associates from pure [AAA] and 
[BBB] -associates. For example, instead of g[AAB]> tne Gibbs free 
energy Ag[ A AB] of the following reaction can be used, 



2 : [AAA]+ 1 -[BBB]^[AAB], 



(11) 



3 L J 3' 
so that 

2 1 

A? [AAB] - 8[AAB] - ^8[AAA] ~ "^8 [BBB] ■ (12) 

The other Gibbs free energies Ag^j are defined similarly. 

Note that Eqs. (5) define «[aaa] and ri[ BBB ] as linear functions 
of the mole numbers of the solution components n A and n B and 
of the mole numbers of the mixed associate types with 
[ijk] ^ [AAA],\BBB]. The mole numbers of the latter associates, 
which consist of both A and B particles, are the internal variables 
of the model. These are determined by minimising the Gibbs 
free energy of the solution at constant «a and n B , subject to the 
mass balance constrains of Eqs. (5). The equilibrium values of 
ti[ij k ], with [ijk] ^ [AAA], [BBB], are found from 



dG 



dn 



= o, 



(13) 



where the derivative in is calculated for fixed ua, n B and 
fixed five variables with [i'j'kf] ^ [ijk] , [AAA] , [BBB] . By 

the chain rule, Eq. (13) reads 



dG dri[AAA] | dG dn\ BBB] ^ dG 



3" [AAA] Mm] 



^ n [BBB] Mm] 



Mm 



o. 



(14) 



Therefore, combining the last equation with Eqs. (5) and (9) 

gives 

&g[AAB] 



n m 



(n[AAA]) am/ \n[ BBB] fm/' 



= exp - 



RT 



(15) 



Here, a^j and stand for the numbers of A and B particles 
in the [ijk] -associate, respectively. For example, OC [aab] = 2 and 
P[aab] = !• 

Recall that Eqs. (5) define «[aaa] as a linear function of ha 
and rt[yjt]'s, and ri[ BBB ] as a linear function of n B and n^'s, 
where [ijk] ^ [AAA] , [BBB] . The chemical potential /ja of the 
solution component A is calculated as follows 



Ha = 



dG_ 

dn A 



dG dn 



ijk] 



[ijk] Mm] 



dn A 



+ 



+ 



dG , 


( dn[AAA] 




dri[ AAA ] 


Mm] 


dn[AAA] 


y dn A 


•"I 

[ijk] 


Mm] 


dtiA 


dG , 


( dri[ BBB ] 




dn[ BBB ] 


Mm] 


Mbbb] 


I dn A 


^E 

m 


Mm] 


dn A 



(16) 



dG dn 



[AAA] 



dn 



+ 



[AAA] dn A 
dG dn[ BBB] 



E 

[ijk] 



dG dn 



[AAA] 



^n [aaa] Mm] 



dG 



dn 



dn 



[ijk] 



dn 



[ijk] 



Mm] 

dn A 



1 cc\ 



Here, the sums are taken over i,j,k —A,B such that [ijk] ^ 
[AAA].[BBB]. Using Eqs. (9) and (14), one verifies that 

dm. 



dG 



l [AAA] 



dn 



[AAA] 



8 A 



(17) 



The chemical potential /jg of the solution component B is cal- 
culated in a similar way. 

Note that [AAB], [ABA] and [BAA] associates differ in the 
spatial arrangement of the constituent particles, though they 
have the same "chemical composition" A2B. The same distinc- 
tion holds for [BBA], [BAB] and [ABB] associates of common 
composition AZ?2- Let us assume, for a moment, that only one 
spatial arrangement of associates is allowed for each composi- 
tion, while the other arrangements are prohibited energetically, 
for example, g [BAA] = g [ABA ] = +°° and g [BAB ] = 8 [abb] = +°°- 
In this case, the proposed formalism reduces to the associate 
species model [3], in which associate species are composed of 
three particles. 

Unlike previous modifications of the associate model, we 
take into account all possible spatial arrangements of particles 
in an associate. If two or more associates of different types are 
spatially symmetric to each other, then their molar Gibbs free 
energies are considered equal. For example, 2-particle asso- 
ciates [AB] and [BA] are symmetric and therefore are endowed 
with equal Gibbs energies g[AB] — 8[ba]- 

In more complex cases, however, Gibbs free energy levels are 
ascribed to associates depending on the spatial arrangement of 
particles in them, so that both energy splits and multiple levels 
may occur for associates of common chemical composition. For 
example, if particle sites in 3-particle associates are arranged 
linearly as shown in Fig. 2(a), then g [AAB ] = 8[baa] while g [ABA] 
can be different. Alternatively, if the particle sites are arranged 
as in Fig. 2(b), the three associate types are all symmetric to 
each other, and hence, g [AAB ] = g[ ABA ] = 8[baa[- 

Note that no assumptions on spatial arrangements of particle 
sites within an associate have been made so far in the frame- 
work of the proposed model. In general, it is impossible to de- 
termine in advance the number of energy levels and their mul- 
tiplicities for the associates of a particular composition. Fur- 
thermore, such associates may have the same Gibbs free en- 
ergy of formation, even if they are not spatially symmetric. 
As a reasonable initial approximation, all associates of a given 
composition can be endowed with the same Gibbs energy of 
formation. This assumption can be refined subsequently in the 
process of thermodynamic model optimisation for real chemi- 
cal systems, if use of several energy levels appears to provide 
a better fit to experimental data. 

For the rest of this section and also in Sections 5 and 6, 
we assume, for simplicity, that the molar Gibbs free energy 
of an associate type is completely specified by its chemical 
composition, 



A 8[AAB] = A 8[ABA] = A 8[BAA] =4g 2 .l, 
A 8[BBA] = A 8[BAB] = A 8[ABB] = 4? 1,2- 



(18) 



Here, the subscript "2,1" signifies that the associate consists 
of two A -particles and one Z?-particle, with "1,2" and similar 




(a) Particle sites are arranged along a (horizontal) line. 
The associates [BAA] and [AAB] are mirror reflections 
of each other about axis a; they are not symmetric, 
however, to [ABA]. 




(b) Particle sites are arranged at the vertices of an 
equilateral triangle. All associates are symmetric to 
each other by rotation about axis b, perpendicular to 
the plane of the triangle. 

Fig. 2. Schematic views of different arrangements of particle sites. 

indices understood appropriately. More general case is consid- 
ered in Section 7. 

Using Eqs. (15) and (18), 



(19) 



n [AAA[ 


= "3,0, 


n [AAB] 


= n ]ABA] 


n [BBA[ 


= n [BAB[ 


n [BBB] 


= "0,3 ■ 



In terms of molar fractions, Eqs. (19) read 

X [AAA[ =*3,0) 

X [AAB] = X ]ABA] ~ X ]BAA] = x 2,l/3, 
X [BBA[ = X [BAB[ = X [ABB] = *l,2/3, 
X [BBB] = X Q3 ■ 

Substitution of Eqs. (19) and (20) into Eq. (9) gives 



(20) 



1 ~7f\ 



G = «3,0g3,0+«0,3g0,3 + «2,lg2,l + «l,2gl,2 

+ /?r (n3,olnx3,o + no,3lnxo3) 



(21) 



+ RT (n 2 ,iln(^| i )+n 



1.2 



In 



*1,2 



The mass balance constraints of Eqs. (5) reduce to 

nA = 3n3,o + 2n2,i+«i,2, 
n B = 3n ,3 + 2n i,2 + n2,i • 



(22) 



The equilibrium values of the associate mole numbers are de- 
termined by 



«2.1 



2/3 1/3 
n 3,0 n 0,3 
«1.2 



1/3 2/3 
n 3,0 n 0,3 



= 3 exp ( — 
= 3 exp 



Ag 



2,1 



RT 

Agl,2 



(23) 



5. Case Ag2,i = and Ag\2 = 

Recalling Eq. (18), consider the situation where there are 
no Gibbs free energy changes of forming different associates, 
that is, Ag2,i = and Agi,2 = 0. In this case, the A and B 
particle species mix ideally. The composition of a randomly 
selected triplet of particles follows the binomial distribution, 
well-known in probability theory; see, for example, Ref. [12]. 
More precisely, the probabilities of choosing an associate of 
particular compositions, or, equivalently, the molar fractions of 
appropriate associates, are 



x ii) = X A' 

x 2 ,\ = 3x\x B , 
xi,2 = 3x^x|, 

X ,3 =x\. 



(24) 



Here, xa and x B are the molar fractions of A and B particles. 
Since the total mole number of associates n tot is equal to (tia + 
n B )/3, then the mole numbers of associates are given by 



«3,0 
«2,1 
«1,2 
«0,3 



3{n A + n B ) 2 ' 
n 2 A n B 

(ha + n B ) 2 ' 
n A n\ 



(25) 



2' 



(ha + n B ) 
3(n A + n B ) 2 



By direct inspection, the mole numbers given by Eqs. (25) sat- 
isfy Eqs. (22) and Eqs. (23) with Ag2,i = and Agi,2 = 0. Since 
the Gibbs free energy G given by Eq. (21) is a strictly convex 
function of the mole numbers of associates, then Eq. (23) has 
no other solution satisfying the mass balance. 

Furthermore, if Ag2, l = and Ag \ 2 — 0, then the molar Gibbs 
free energies of associates are given by 



£3,0 =3gA, 

g2,l = 2g A +g B , 

(26) 

gl,2 =gA + 2g B , 

go,i = 3g B . 

Finally, substitution of Eqs. (25) and (26) into Eq. (21) yields 

G = n A gA+n B gs+RT (n A \nx A + n B \nx B ) . (27) 

Thus, the proposed model correctly reduces to the ideal solu- 
tion model in the case where there are no interactions between 
mixing particles. 

6. Dilute solutions 

Pelton et al. [10] pointed out another interesting feature of 
the associate model that occurs in dilute solutions. They consid- 
ered the associated solution of monoparticles A and B and as- 
sociates A2B. The highly ordered solution rich in component B 
consists primarily of B-monoparticles and A2#-associates. Ac- 
cording to the authors, the chemical activity a B of the compo- 
nent B behaves asymptotically as (1 — xa/2) rather than (1 — xa) 
for small xa- Thus, lim (da B /dxA) = —1/2. Note, however, 

that this behaviour of a B is observed only in the limiting case 
Ag2 1 = — °°- For any finite value, lim {da B jdxA) — — 1. We 

x A ^0 

present the proof of this result for the model proposed in this 
paper. Using the presented technique, similar result can be es- 
tablished for the example considered by Pelton et al. [10]. 

1 /3 

As obtained in Section 4, a B = % 3 . The required derivative 
is then calculated as follows. First, 



da B 

dXA 



1 -2/3 ^0.3 
3*0,3 



1 _2/3 / dx A 



dxA 3 0,3 \dxQ^ 



-1 



(28) 



Secondly, assuming Ag2,i an d Ag\ ,2 finite and denoting the right 
hand sides of Eqs. (??) by 



£2,1 = 3 exp 
we obtain 



Ag2,l 

RT 



Ei, 2 = 3 exp 



2/3 1/3 



As 



1.2 



RT 



X2,l =£2,1X 3 X 3 , 



1/3 2/3 

Xl,2 ^1,2X3 qXq 3 . 



(29) 



(30) 



Substitution of Eqs. (30) into the mass balance constraints of 
Eqs. (22), with the latter written in terms of molar fractions, 
gives 

2/3 1/3 1/3 2/3 

1 =*3,0 + *0,3 + £2,1*3,0 *0,3 + e l,2*3,0 x 0,3 ' ( 31 ) 
1 — Xa — X ,3 + 3 e 2,l x 3.0 x 0,3 + 3 £ 1,2 X 3,0 X 0.3 - (iZ > 

Now, Eq. (31) implicitly defines xo,3 as a function of X3 0, while 
Eq. (32) defines xa as a function of X3,o and ^0,3. Differentiating 
Eq. (31) gives 



dxn-i 



fei,2f- 



\z2.\t 2 



l + 5ei,2f" 1 + le2,if" 2 ' 



(33) 



1 71 



where 



X3£ 
*0,3 



1/3 



(34) 



Differentiation of Eq. (32) with respect to xo,3 an d substitution 
of Eq. (34) into the resultant expression yields 



dXA 

dx Q 3 



1 £i it 

9 *> z 



1 



7E2,ir 



(35) 



2 (ei, 2 ?~ 2 + e 2 ,i?~ 1 ) (l + fei,2f + |£2gf 2 ) 
"9 (i + i ei ^-2 + 2 e2 ^-i) 

, 4 1 2 

= -l-^e li2 f--e 2 ar 2 

2 (ei : 2 + £2,i0(l + fei,2? + ^£2,if 2 ) 

9 (f 2 + |e li2 + |e 2: if) 

Note that xoj — ► 1, X3_o — > and f — > as — > 0. Therefore, 
substituting Eq. (35) into Eq. (28) and taking the limit as xa — > 
0, one verifies that lim (dag/dxA) = —1 for any finite values 

* A ->0 

of Ag 2: i and Agi j2 . 

Now consider the limiting case where Ag 2 ,i — > — °°, while 
Agi2 remains finite. Recalling Eqs. (29) and taking the limit in 
Eq. (35) give 



lim 



dXA 



1 



1 



62,1" 



= El 2 f 1 . 

•dxQT, 3 6' 6 



Finally, substituting Eq. (36) into Eq. (28) and taking the limit 
as xa — > 0, we obtain lim {dag/dxA) = —1/2 as Ag 2 i — > — oo. 

x A ^0 



7. Model equations for arbitrary number of solution 
components and arbitrary size of associates 

Consider an r-component solution A \ — . . . — A r and assume 
that the components and the solution itself all consist of m- 
particle associates. There are r m distinguishable types of asso- 
ciates. Some of these types have the same compositions. Now 
consider an m-particle associate that consists of k\ , . . . , k r par- 
ticles of types A\,...,A r , respectively. Thus, the composition 
of the associate is specified by the r-tuple of nonnegative inte- 
gers (fei , . . . , k r ) satisfying k\ + . . . + k r = m. Omitting the de- 
pendence on r, the set of such tuplets, which represent all pos- 
sible compositions of m-particle associates in the r-component 
solution, is denoted by & m . Its cardinality, that is, the number 
Ncomp of different compositions is computed as 

(m + r-l)! 



m\(r-l)\ ' (3?) 

Consider a pure solution component A, . Its m-particle associates 
all have the composition 

o, = ((),..., 0,m, (),..., 0). (38) 

i— 1 r— i 

The compositions G\,...,O r and corresponding associate 
types are referred to as pure. The complementary set 
&m\{Gi, - ■ ■ ,G r } of mixed compositions, containing two or 
more different particle species, is written briefly as &° m . Thus, 
&^ is constituted by those r-tuples (ki,.. .,k r ) from & m with 
at least two nonzero entries. 



In one mole of the solution component A,-, there are 1 /m 
moles of m-particle associates of composition a,. Thus, their 
molar Gibbs free energy g aj is 

goi = mgt, (39) 

where gi is the molar Gibbs free energy of the solution com- 
ponent A;. 

In general, the total number N kl kr of distinguishable types 

of associates which have the same composition ki 7 ...,k r is 
described by the multinomial coefficient 

m! 

N kl _ kr = . (40) 

1=1 

We assume that all distinguishable types of the associates of 
composition (ki,...,k r ) are endowed with J kl different val- 
ues of the molar Gibbs free energy. Let d k J kr be the number 

of the associate types at the jth energy level g k j ^, that is, 
the multiplicity of the level, so that 



\,...* r 



E < kr =N ku .., kr . 

7=1 



(41) 



Qg-j All the d k j kr associate types have the same molar Gibbs free 



energy of formation Ag[ j j kr according to the reaction 

^ (A , ) m + . . . + (A,.),„ +± ( (A , ) kl . . . (A r ) kr ) W . (42) 
m m 

More precisely, the Gibbs free energy of the reaction is defined 

by 

a4!.^=4U-^-----^- («) 

The molar Gibbs free energies g^J k of mixed associates of 
composition [k\ , . . . , k r ) € &° m , or alternatively, the correspond- 
ing formation energies Ag k j j ^ from Eq. (43) are adjustable 
parameters of the model. 

Now let r$ , denote the mole number of the associates 

k\,...,k r 

of composition (k\, . . . ,k r ) at the y'th energy level. The mass 
balance constraints then read 



■V...,t r 



(44) 



»«= E k i E % & i = 1 > 

{h,...,k r )e& m i=i 
where n, is the mole number of A, particles in the solution. One 
verifies that 

r J k l ,-,k r 

£n,=m £ £ n ki..,k r = mntot - (45) 
(=1 (ki,...,kr)e&„ j=i 

Note that the total number of associates n tot is independent of 
composition of the solution since they are assumed to be of 
equal size. The Gibbs free energy of the solution is computed 

as 



h\ k r 



E E 4! : J, *, 

(k u ...,k r )e& m ;'=i 



Jk\ ,...,kr s ■ ' 

+ RT E E < k, ln 



(46) 



(k 1 ,..., kr )ee m j=i 



d 



[j] 

k\,...,k r 



1 70 



where 



%,...,k r 



1 k 1 ,...,k r 

«tot 



(47) 



is the molar fraction of an appropriate associate type. 

Eqs. (44) define n 0| as a linear function of the mole number 
m of the solution component A, and of the mole numbers of the 
mixed associates n^j ^, where {k\,...,k r ) G &° n . The equi- 
librium values of n^j ^, which are internal variables of the 
model, are determined by minimising the Gibbs free energy of 
the solution at constant m,...,n r , subject to the mass balance 
constraints of Eqs. (44). The minimum is found by setting 



dG 



dnf , 

ki,...,k r , 







(48) 



»l,...,n r 

for all (ki,...,k r ) G &° m and for all [j] = l,...,J kl kr . Recall- 
ing that n Gj are functions of n^j ^, where (ki,...,k r ) G &° n , 
and using the chain rule, 



dG ^ dG dn a . 



= o. 



Substitution of Eqs. (44) and (46) into Eq. (49) gives 

J;1 



hr^k - J j] exn I - 
-%,...,k r ex P 



i=i 



°h,..., k r 

RT 



(49) 



(50) 



Using Eqs. (44), one verifies that n 0i is a function of n, and of 

"if k r ' wnere (^i> • • • >kr) € 6° r The chemical potential /ji of 
the solution component A i can now be calculated as follows: 

dG' 



A<i = 



dn 



n 2 ,...,n r 

4, * _dG dn l l 
L . A- a J/ 



k\,...,k r 



(k u ...,k r )e&° m .7=1 3njt,,...^ 



3«i 



+ E 



3G / dn 0i t y ''Y' <K,. H ... 



^dn a , \ dm {h _ kr)e&m p x dn yj_ kr 



dG dn 



CM 



dn a , dm 



^ ^ dG dn Gl 



r<K L 



dG dn 




(ki,...,k r )e6° m j=l \dn 



h,...,k r 



dn Gl dm 



(51) 

Here, Eq. (49) have been used. Using Eqs. (39) and (44), one 
verifies that 

1 



A*i = gi + —RTlnx 0l 



(52) 



The chemical potentials of the other components are calculated 
in a similar way. 

Now assume that all the associate types of composition 
(k\ , . . . , k r ) have the same Gibbs free energy of formation 
for any (k u ...,k r ) G & m , so that Jk u ...,k r = 1 and d\ u K = 



m\/Y\kj\. This assumption, which can be subsequently re- 

fined in thermodynamic model optimisation of real chemical 
systems, allows substantially reduce the number of adjustable 
parameters of the model. Since the Gibbs free energies of pure 
associates are fixed, the number of adjustable parameters is 
equal to (N comp — r), where N comp is given by Eq. (37). Then, 
omitting the superscript [j], Eq. (46) reduces to 



E n *i,.. 

(k l ,...,k r )£<S m 



8k, 



RT\n 



Xk u 



1=1 



\ \ J 

(53) 

Similarly to Section 5, consider the case when all Gibbs free 
energies Ag^ ..^ of formation of associates are equal to zero, 
so that particles of the solution components are mixed ideally. 
In this case, the composition of m randomly selected particles 
follows the multinomial distribution; see Ref. [12] for more 
details. The molar fraction of the associates with composition 
(k\ , . . . ,k r ) G & m is expressed as 



i=i k '- 



and hence, their mole number is 

(m-l)l 



n k k, 



L li- 
fe i 



1=1 



(54) 



(55) 



Here, n, and x, are the mole number and the molar fraction of 
the solution component A,. One verifies that the mole numbers 
given by Eq. (55) describe the solution of Eq. (50) subject to the 
mass balance constraints of Eqs. (44) in the case Agk h ....k r = 0. 
The uniqueness of the solution is ensured by the strict convexity 
of the Gibbs free energy G given by Eq. (53) in the variables 
"I. <.• 

From Eqs. (39) and (43), the molar Gibbs free energy of the 
associates of the composition (ki,.. .,k r ) is given by 



gk u ...,k r = Y< ki 8> ' 



(56) 



where, gi is the molar Gibbs free energy of the solution com- 
ponent A,. Substitution of Eqs. (54)-(56) into Eq. (53) and a 
straightforward, though lengthy, verification shows that the pro- 
posed model correctly reduces to the r-component ideal solu- 
tion model in this case. That is, 



G=Y d n i (g i +RT\nx i ) 



(57) 



i=i 



From this reduction, it is immediately follows that the model 
with associates of size r can be reproduced by the models with 
associates of size 2r, 3r and so on. For example, consider the 
model with associates of size 2r. Any 2r-associate is a combi- 
nation of two r-associate. Now assume that there is no Gibbs 
free energy change on forming any 2r-associate from corre- 
sponding two r-associates. As demonstrated above, the model 
with 2r-associate reduces to the ideal mixture of r-associates. 
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Thus, the model with 2r-associates is more general and include 
the model with /--associates as a particular case. 

8. Effective adjustable parameters of the model 



Gibbs free energies of associates, the optimal values of the ef- 
fective adjustable parameters k r can be determined by the 
trial-and-error procedure which is conventionally used in ther- 
modynamic model optimisation of real chemical systems. 



The adjustable parameters of the model related to asso- 
ciates of the composition (k\,.. .,k r ) are Ag^j ^, where [j] = 
!)• • ,Jki,...,k r - One should also take into account the multiplic- 
ity of the energy levels d^j k . In general, the number of ad- 
justable parameters of the model increases exponentially with 
the increase in size of associates m. However, the number of 
adjustable parameters can be substantially reduced without loss 
of generality of the model as described below. 

Consider associates with the composition (k\,...,k r ). In gen- 
eral, Jk\,...,k r energy levels are possible for these associates. The 

mole number n/ Cl k r of all associates with the composition 

(k\,.. .,k r ) is given by 

\,...,kr 

(58) 

7=1 

The mass balance constraints of Eqs. (44) take the form 

m= E k '"t: i=l,...,r. (59) 

(ki,...,k r )ee m 

Summing up over [j] in Eq. (50), one verifies that the equilib- 
rium value of riky k r is calculated as 

.,k r 



n h ,...,k r = E n kl..,k r 



= z 



where 



nK) fei/m 
i=i 



h x ,....k r 



k\ ,...,k r 



(60) 



7=1 



RT 



(61) 



The function Zk y .....k r defined by Eq. (61) plays a role of the 
partition function which describes the distribution of associates 
of the composition (k\ , . . . ,k r ) over energy levels. That is, 



J>1 



»£, ky 



Xk u 



■ exp - - 



Zki,...,k r 



RT 



(62) 



In fact, Zk x k r is a single effective adjustable parameter related 

to associates of the composition (k\ , . . .,k r ). The other thermo- 
dynamic parameters of the model can be expressed in terms of 

Zk x k r , where k\,... 7 k r e &,„, and their derivatives. Indeed, 

using Eqs. (39), (43), (59) and (62), one verifies that the Gibbs 
free energy of the solution given by Eq.(46) takes the form 



G = E" 1 *' 

;=i 



RT 



E b *i. 

..,k r )e6 m 



In 



Xk u ...,k r 

Zk, 



-A 



(63) 



According to Eq. (61), Zk h ... t k r varies from zero to infinity. 
When the associates of the composition {k\,.. .,k r ) are prohib- 
ited energetically, that is (Ag^j ^ — > +°°) for all [j], Z^ k r 

approaches zero. If the associates of the composition (k\ , . . . , k r ) 
is highly preferable, k r approaches infinity. Similar to the 



9. Excess Gibbs energy terms 

Similarly to the modified associate species model [4], regular 
or, in general, polynomial, excess Gibbs free energy terms can 
be included into the proposed formalism. These terms take 
into account interactions between associates. A probabilistic 
interpretation of the polynomial excess Gibbs free energy terms 
is presented in Ref. [13]. In fact, this interpretation provides a 
theoretical justification for such terms. Treating the associates 
as particles, the results of Ref. [13] are applicable to the model 
presented in this paper. 

Note that it is desirable to use the excess terms only for "fine 
tuning" of the model, while the main adjustable parameters are 
the molar Gibbs free energies of associates. The following two 
conditions on the absolute value of the interaction parameters 
should be satisfied; see, for example, the monograph by Pri- 
gogine and Defay [8] for more details. 

(i) The absolute values of the interaction parameters should 
be small compared with the molar Gibbs free energies of 
the associates. As pointed out by Prigogine and Defay [8], 
if the interaction between associates, say "CI" and "C2", 
is sufficiently strong to alter the vibrational and rotational 
states of the associates, then the associate "C1C2" is 
included into the set of associates by the definition of 
the classical associate model. In the framework of the 
proposed formalism, one should consider the associates 
of larger size. 

(ii) The values should also be small in comparison with RT. 
Otherwise, the assumption of ideal mixing of associates 
is less justified. Again, larger associates should be taken 
into account. 

Note, however, that the second condition is sometimes re- 
laxed in order to fit the experimental data available for real 
solutions. This is the case, for example, for the solutions with 
immiscibility, which is the result of relatively weak, compared 
with the Gibbs free energies of associates, repulsive interac- 
tions between the associates. 



10. Discussion on applicability of the model 

The suggested modified associate formalism belongs to 
associate-type models. As a result, the range of applicability 
of the formalism is at least the same as that of the classical 
associate model or the associate species model. There are, 
however, some distinctions in using the suggested formalism 
and the previous modifications of the associate model. These 
distinctions are discussed below. Since the present study is 
intended as a theoretical introduction to the modified asso- 
ciate formalism, the results on its application to real chemical 
systems will be reported elsewhere. 
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In contrast to the previous modifications, where associates 
of arbitrary compositions can be included into the set of asso- 
ciates, the compositions of associates are defined by their size 
in the framework of the modified associate formalism. There- 
fore, a modeller should pay special attention to selection of the 
size of associates. If experimental information about composi- 
tions of associates that present in the solution phase is avail- 
able, this information should clearly be taken into account. The 
size of associates should be large enough to incorporate those 
compositions. 

When the suggested formalism is applied for thermodynamic 
description of multicomponent solution phases, the size of as- 
sociate should be large enough to incorporate the associates of 
the compositions, which coincide with those of maximum or- 
dering in all binding binary systems. For example, consider the 
ternary system A — B — C and assume that in the binary system 
B — C the composition of maximum ordering is that of the as- 
sociate B2C, while in the systems A—B and A — C maximum 
ordering occurs at the compositions of the associates AB and 
AC, respectively. To incorporate the required composition, the 
size of associates should be divisible by 2 and by 3. Therefore, 
the size of associates for the ternary system A— B — C should 
be at least 6. 

It is desirable to use 6-particle associates for thermodynamic 
model optimisation for all the binary systems. There is, how- 
ever, no need to use 6-particle associates from the beginning. 
Let us assume that the systems A-B and A-C are initially op- 
timised with 2-particle associate, while the system B-C is op- 
timised with 3-particle associate. As discussed in the end of 
Section 7, the model with associates of size r can be repro- 
duced by the models with associates of size 2r, 3r and so on. 
Using this property of the proposed formalism, the descriptions 
of the binary systems with 2- and 3-particle associates can be 
replaced by the equivalent descriptions with 6-particle asso- 
ciates in a straightforward way. Then, the binary systems can 
be combined into the ternary one. 

As a demonstrational example, consider the system A — B. 
When the system is optimised with 2-particle associates the 
Gibbs free energies of the associates g[ AA ], g[BB]i 8{ab] and 8[ba] 
are known. Due to spatial symmetry, g^g\ = 8[ba\- In the equiv- 
alent description of the system A—B, 6-particle associates are 
formed from three 2-particle associates with no Gibbs free en- 
ergy changes on such formations. The associate of the composi- 
tion A& is formed from three [AA] -associates with gA 6 — 3g[AA]- 
The associate of the composition A 5 B is formed ether from two 
[AA] -associates and [AB] -associates or from two [AA] -associates 
and [BA] -associates. Then, recalling that g^AB] = 8[ba]> A5B- 
associates have only one energy level gA 5 B — ^g[AA] + 8[ab] with 
multiplicity 6. The A4B2-associate can be formed from two 
[AA] -associates and one [BB] -associate with the energy level 

8aIb 2 = 2 8[aa] +8[bb] and the multiplicity d [ X = 3. Alterna- 
tively, A4Z?2-associate can be form from one [AA] -associate and 
either two [AB] -associates or two [BA] -associates or one [AB]- 

associate and one [BA] -associate. In this case, the energy level 

[2] [2] 
is 8a 4 b 2 = 8[aa] +2S[as], and its multiplicity is d\\ Bi = 12. In a 

similar way, one verifies that A3B3-associates have two energy 



levels: g^ = g [AA ] +g[ B B] +8[ab\ with the multiplicity d£ Bj = 

[2] [2] 
12 and g A ^ 3 = 3g[ A B] with the multiplicity d\^ = 8. A 2 B 4 - 

associate also have two energy levels: = g\AA] + 2 8[bb\ 

with the multiplicity = 3 and g®^ = g [BB] + 2g [AB] with 

[2] 

the multiplicity d A % =12, while AB5 -associates have only one 
energy level g A B 5 = ^g\BB] +8[ab] with the multiplicity equals 
to 6. Bg-associates also have only one energy level gs 6 — 3g^ B B] 
with the multiplicity equals to 1 . The binary systems A — C and 
B — C are treated similarly. 

Previous modifications of the associate model allow more 
flexibility in fitting experimental data compared with the sug- 
gested approach, since stoichiometry of associates can be ar- 
bitrarily altered to improve the fit. This additional flexibility 
can be helpful for modelling a binary system with high degree 
of ordering, where the described paradox is of relatively low 
practical importance. However, when such a binary system is 
combined with other binaries to form multicomponent model, 
arbitrary selection of associate stoichiometry can result in the 
negative consequences of the entropy paradox. 

Similar to the previous modifications, a miscibility gap can 
not be reproduced without excess Gibbs free energy terms in the 
proposed formalism. However, as demonstrated by Besmann et 
el. [4] for the associate species model, immiscibility, which is 
the result of repulsive interactions between associate (or asso- 
ciate species), can be accurately represented by the associate- 
type models with the polynomial excess Gibbs free energy. 
Treating associates (or associate species) as particles and using 
the results of Ref. [13], one verifies that the coefficients of the 
polynomial excess Gibbs free energy explicitly relate to ener- 
gies of interactions between associates. It is important to note 
that, in the case of the modified quasichemical model [10], im- 
miscibility is described by empirical polynomial expansions of 
the Gibbs free energy of the quasichemical reaction. In contrast 
to associate-type models, physical meaning of the coefficients 
of such expansions is not clear. 

Another distinction of the proposed formalism from the pre- 
vious modifications is related to multicomponent associates, 
that is, the associates which consist of particles of three or more 
types. In the previous modifications, no multicomponent asso- 
ciates are initially considered. This implies that all multicompo- 
nent associates are assumed to have infinite positive Gibbs free 
energy. Multicomponent associates are usually included into the 
model, when such associates are required to fit available exper- 
imental data in multicomponent systems. In contrast, the com- 
positions of multicomponent associates are prescribed by the 
selected size of associate in the proposed formalism. It would 
be beneficial to develop a method for estimating the Gibbs free 
energies of the multicomponent associates from those of bi- 
nary associates. The development of such a method, however, 
is the topic for separate study. As an initial approximation, the 
Gibbs free energies of formation of the multicomponent asso- 
ciates from pure associates can be set to zero. If experimental 
data are available in multicomponent systems, the Gibbs free 
energies of the multicomponent associates can be adjusted in a 
conventional trial-and-error procedure. 
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A possibility to select compositions of multicomponent as- 
sociates arbitrarily, which exists in the previous modifications, 
could provide more flexibility in data fitting compared with the 
fixed set of compositions. It could be the case, that arbitrarily 
selected additional multicomponent associates provide better 
fitting to particular experimental data. One should recognise, 
however, that the price for the possibility to select the composi- 
tions of associates arbitrarily is the entropy paradox, which un- 
dermines the fundamental predictive capabilities of the model. 
In our opinion, the necessity of additional multicomponent as- 
sociates for fitting the available experimental data indicates that 
the size of associates should be increased to incorporate the 
required compositions. 

11. Conclusion 

In this paper, the Modified Associate Formalism has been 
proposed for thermodynamic modelling of solutions. The pre- 
sented approach is free from the entropy paradox and correctly 
reduces to the ideal solution model, where there are no Gibbs 
free energy changes on forming associates of different types. 

Asymptotic behaviour of chemical activities of solution com- 
ponents in binary dilute solutions has been investigated. It has 
been demonstrated that the derivative of the chemical activity 
of a solution component in its molar fraction at terminal com- 
position has the expected value for any finite value of the Gibbs 
free energies of formation of associates. 
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